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Abstract 

We demonstrate that chaos can be controlled using a multiplicative exponential feed- 
back control. All three types of unstable orbits - unstable fixed points, limit cycles and 
chaotic trajectories can be stabilized using this control. The control is effective both for 
maps and flows. The control is significant, particularly for systems with several degrees 
of freedom, as knowledge of only one variable on the desired unstable orbit is sufficient 
to settle the system on to that orbit. We find, that in all the cases studied, the transient 
time is a decreasing function of the stiffness of control. But increasing the stiffness beyond 
an optimum value can increase the transient time. The control can also be used to create 
suitable new stable attractors in a map, which did not exist in the original system. 
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1 INTRODUCTION 

The problem of controlling chaos has recently received much attention [1-8] . A chaotic 
system in general cannot be made to converge to a freely evolving desired trajectory, 
whether periodic or chaotic, because of the inherent unpredictibility of the system. The 
control of chaos in this context consists of forcing the system to evolve along a desired 
trajectory. Pecora and Carroll [1-3] demonstrated that two identical chaotic systems driven 
by a common signal display asymptotic convergence of their trajectories even though they 
may have started from very different initial conditions provided the Lyapunov exponents 
of the driven system are all negative. Such behaviour was demonstrated numerically and 
also confirmed experimentally. 

Ott, Grebogi and Yorke (OGY) [4,5] succeeded in forcing a chaotic system on to one of 
its own unstable periodic orbits by making a set of small time-dependent perturbations on 
the system parameters. They demonstrated their method numerically by controlling the 
Henon map. They also pointed out that as an infinite number of unstable periodic orbits 
are embedded in a chaotic attractor, the stabilization of these unstable periodic orbits can 
lead to an enhancement of the system's performance as well as to the adaptibility of the 
system to varying performance requirements since the system behaviour can be changed 
by perturbing the system to stabilize different unstable orbits. In the absence of chaos, 
separate systems would be required for each of the different responses. The effectiveness 
of the OGY method of control has been confirmed experimentally in several systems [6-8]. 

In this paper we have succeeded in implementing a novel method for controlling chaos. 
We have used a multiplicative exponential feedback control on a parameter of the system, 
with the argument of the exponential being proportional to the feedback response of the 
system, i.e. the difference between the desired value and the actual value of one of the 
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suitably chosen variables of the system. Our method is a combination of variable and 
parameteric control and it can stabilize all three types of unstable orbits - unstable fixed 
points, limit cycles and chaotic trajectories. It is found to work effectively both for maps 
and flows. 

Consider a general N-dimensional dynamical system 

t = tQt-i,-t) (1) 

where it = {Xi,X2, , -^jv) are variables and = (//i, fi2, I^k) are parameters whose 

values determine the nature of the dynamics. The stabilization of a desired unstable 
attractor or a chaotic trajectory is possible by multiplying a suitably chosen parameter, say 
in (1) by an exponential feedback control involving only one suitably chosen variable, 
say Xi, with the form of control being given by 

exp[e{Xi - Xt)]. (2) 

Here Xi is the actual value of one of the variables of the system after applying the control, 
Xf is the desired value of that variable and e is the "stiffness" of the control which can 
take both positive and negative values. 

The dynamics of the modulated system in the presence of control is given by 

It = 'P{Jt; III, ii2, iirexp[e{Xi - Xf )], .., /ik] t). (3) 

Note that the control becomes passive once the desired goal X^ is achieved. If fluctuations 
drive the system off the desired goal, the control reactivates. 

The control works for those combinations of controlling parameters and variables of 
the system for which the largest real part of the Lyapunov exponents of the modulated 
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system, represented by eq. (3), is negative. The feedback function in the expression of 
the control involves only one suitably chosen variable Xi to convert the desired repellor 
- whether a fixed point, a limit cycle or a chaotic orbit, into an attractor. This shows 
that the knowledge of only one variable on the desired unstable orbit is sufficient to settle 
the system on to that orbit. This makes the control particularly useful for systems with 
several degrees of freedom. For example, a desired unstable fixed point of the Lorenz 
system [9] can be stabilized using control with feedback depending on any one of the X, 
Y or Z variables whereas an unstable limit cycle or a chaotic trajectory of the Lorenz 
system can be converted into an attractor using feedback depending on the Z variable 
only. For stabilizing most repellors, it is sufficient to multiply one chosen parameter by 
the exponential control; but sometimes, as in the case of one of the two fixed points of the 
Henon map [10], both the parameters a and b of the system have to be multiplied by the 
exponential control. Although in all the cases we have studied, only one variable on the 
desired orbit is sufficient to make all the Lyapunov exponents of the system negative, it 
may not be possible in some other systems where more than one variable may be required 
for the stabilization of unstable orbits. But the fact remains that the control uses only 
a subset of the variables and the parameters for controlling chaos. We have tested our 
control for stabilizing different types of orbits for the logistic map, the Henon map and the 
Lorenz system and found it to work effiectively. 

A quantity of obvious interest in the context of controlling chaos is the time required 
for the system to settle on to the desired orbit. This of course depends upon the stiffness 
of the control i.e. e. For a given e we study the time r required for the system to approach 
within a distance uj of the desired orbit starting from some initial point. If uIq is the initial 
distance from the desired orbit, then it is clear that the length of the transient r and u 

4 



are related by a; = u!oexp{XT), where A is the largest real part of the Lyapunov exponents. 
The slope of the plot of r against ln{u>/u>o) is nothing but 1/A. This of course has to 
be negative for convergence. The values of the Lyapunov exponents calculated from such 
points are found to be in good agreement with those obtained either numerically using 
the method given in [11] or analytically wherever this is possible. We have studied the 
transient time r required for settling on to a desired orbit within a given accuracy a; as a 
function of e and found that r in general is a decreasing function of e. But there exists an 
optimum stiffness of control beyond which increase in e increases r. This behaviour of r 
is found to be the same as that of the A with e for any orbit. 

Finally, for discrete maps we have also tried to stabilize orbits which are not the natural 
fixed points (stable or unstable) of the original system. We have found that the control 
succeeds in creating desired new stable attractors which are not the natural attractors of 
the unmodulated system, depending on our requirement and use, using the discrete map. 
However we cannot stabilize any arbitrary orbit. The functional form of the map and the 
control criterion decides which orbits can be forced on to the system. 

There is one drawback. The form of the control is such that for a given system 
the unstable fixed point represented by a null vector cannot be stabilised. The reason is 
obvious. The Jacobian of the modulated system given by (3) evaluated at such a point 
in the presence of the control is the same as that of the unmodulated system (1). So the 
eigenvalues of the system remain unchanged in the presence of the control. Hence such a 
fixed point remains unstable even under control. Such points may, however, get stabilized 
in the process of stabilizing other fixed points, resulting in the coexistence of more than 
one attr actor. 

The organization of this paper is as follows. In section 2 we study exponential control 
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for stabilizing an unstable fixed point. We deal first with 1-D discrete systems, then with 
2-D discrete systems and finally with continuous dynamical systems. In section 3.1 we 
extend the control algorithm so as to stabilize unstable closed orbits for flows. We suggest 
a more effective control for stabilization of higher orbits of a discrete dynamical system in 
section 3.2. In section 4 we discuss the stabilzation of chaotic orbits. We use our control 
for stabilization of arbitrary fixed points for discrete systems in section 5. Finally we 
summarize our results in section 6. 

2 EXPONENTIAL CONTROL FOR STABILIZING AN 
UNSTABLE FIXED POINT 

2.1 1-D Discrete System 

A representative one dimensional map is the logistic map given by 

Xn+l = 4/iX,(l - Xn) = Fill, X,) (4) 

where < X, < 1. The map has two period one fixed points corresponding to any 
value of 11, which are X* = and X^ = 1 — 1/(4//). The fixed point X* is unstable for 
/X > 0.25 while X^ is unstable for both < /x < 0.25 and /x > 0.75. As mentioned in the 
introduction, the control is ineffective in stabilizing the unstable fixed point X*. 

For stabilizing X^, we multiply /x by exp[e{Xn — X^)] in (4) so that the logistic map 
in the presence of control is given by 

Xn+i = 4fiexp[eiXr, - X;)]Xnil - X,) = J^ifi, Xn). (5) 

For a given /j, and for e in the range emin < e < ^maxj the Lyapunov exponent is negative 
(i.e. ^Ix* liss in the interval (—1,1)) and hence X^ becomes a stable fixed point (an 
attractor) of the modulated map. (This is found to be a general feature of the control 
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for maps as well as for flows.) The actual expressions for emin and Smax for (5) are: 

Cmm = 4//(4// - 3)/(4// - 1) and Smax = 4//. 

For a given /j, and u, the transient r is found to be the decreasing function of e, 
however, there is an optimum stiffness control, beyond which increasing e increases the 
transient time. This behaviour is also seen in the variation of the A with e. The optimum 
stiffness of control corresponds to that value of e for which A is minimum. For the logistic 
map, A = — oo for e = 4// (4// — 2)/ (4// — 1) and the system settles to more slowly both 
for larger and smaller e. 

2.2 2-D Discrete System 

We have studied the Henon map [10] as an example of 2-D maps. It is given 

by 

Xn+i ^Yn + l-aXl 

(6) 

where X and Y are variables and a and b are the controlling parameters. The map has 
two critical points {X±,Y±) where 

X± = {2a)-\-l + b ± V(l-^)'+4a) 

y± = bx± 

These two critical points are unstable for a = 1.4 and b=0.3. 

To stabilize 1+), we multiply the parameter b by exp[e{Xn — X+)] so that the 
dynamics in the presence of control is given by 

Xn+l = Yn + 1 — aX^ 

(7) 

Yn+i = bexp[e{Xn - 

The other unstable critical point becomes stable only when both parameters a and 
b are multiplied by exp[e{Xn — X-)], so that the map in the presence of the control for 
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stabilizing {X_,Y_) is given by 

Xn+i = Yn + l- aexp[e{Xn - X-)]Xl 

(8) 

Yn+i = bexp[eiXn-X^)]Xn. 
For a given value of a and b, the behaviour of the transient time r with e for a given 
accuracy u> is similar to that seen in discrete 1-D maps. 

2.3 Continuous Dynamical System 

We have studied the Lorenz system [9] as an example of a continuous dy- 
namical system. This is governed by the equations 

X = a(Y-X) 

Y = -XZ + rX-Y (9) 
Z = XY -bZ 

where X, Y, Z are variables and a, r, b are the controlling parameters. The system has 
three critical points, viz 

X' = 0, y = 0, Z' = 
X" = ^yb{r-l), Y" = ^yb{r - 1), Z" = r - 1 

and 

X'" = -Vb^^l Y'" = -^M^l Z'" = r - 1. 

We choose cr = 10, 6 = 8/3 and r — 60 for which all three critical points are unstable fixed 
points of the Lorenz system. The first critical point {X' ,Y' , Z') cannot be made stable 
using our form of control. This has been discussed in the introduction. The second and 
third critical points can be stabilized by multiplying b by exp[e{Z — Z")] and exp[e{Z — 
Z'")] respectively. It is also possible to stabilize these points using other combinations of 
parameters and variables, provided the control criterion i.e., the largest real part of the 
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Lyapunov exponents, A, is negative, is satisfied. Such combinations, with the corresponding 
values of A, are hsted in table I for e = —0.001, —0.01 and —0.1. Similarly, A can be worked 
out for positive values of e also. The table shows that the r-X, r-Y and b-Y forms of control 
work for suitably chosen negative values of e. On the other hand, it is the r-Z, b-X and b-Z 
forms of control which are effective for suitably chosen positive values of e. For the rest of 
the forms viz. cr-X, cr-Y and cr-Z, the value of A remains positive for all e (the Lyapunov 
exponents with these forms of control are independent of e and are the same as that of 
the unmodulated system) thereby implying that these forms of exponential control cannot 
stabilize the unstable fixed points of the Lorenz system. 

For a fixed u>, the transient time r, using b-Z exponential control, is found to be 
a decreasing function of e. However there exists an optimum value of e beyond which 
increasing e increases r (as in the case of discrete 1-D maps). The Lyapunov exponents 
obtained analytically for the second and third critical points in the presence of control 
have the same value. We plot A vs e in fig. (1). The variation of A with e shows the 
same behaviour as that of r with e. The optimum value of e for which r is minimum 
also corresponds to the most negative value of A and is 0.2 for the given values of the 
parameters. The value of A obtained analytically is —6.8993 for e = 0.2. On plotting r, 
as a function of ln(a;/a;o) for the same value of e, we find that the reciprocal of the slopes 
for the second and third critical po! ints are ~ —6.92, which is in good agreement with 
the value obtained analytically. To compare our results with those of Pecora and Carroll, 
with both X and Y drives, we first consider the control with the parameter r multiplied by 
the exponential feedback function involving the X variable for the same parameter values 
(such a choice of the control is made because Pecora and Carroll's method uses variables 
X and Y as drives). We find that the value of A with this form of control is —4.536 for 
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e = —0.17. The value of A for control involving the parameter r and the variable Y is 
—7.593 for e = —0.05, for the same parameter values. The corresponding values of A in 
the method suggested by Pecora and Carroll for the same parameter values and are —1.83 
for the X drive and —2.85 for the Y drive [12]. Thus our values of A are more negative 
and hence the control asserts itself quicker. This can be seen in fig. (2) for the case of the 
second unstable fixed point of the Lorenz system. This result is not surprising in view of 
the fact that our control has an exponential form. 

3 EXPONENTIAL CONTROL FOR STABILIZING 
UNSTABLE HIGHER PERIOD ORBITS 

3.1 STABILIZATION OF UNSTABLE CLOSED ORBITS 
FOR FLOWS 

The unstable limit cycle of a continuous system can be converted into a stable attrac- 
tor by multiplying one suitably chosen parameter, say /ii, (or may be more parameters 
depending on the system) by the exponential feedback function which depends only on 
one of the variable say A^, i.e. 

fxi iJ,iexp[e{Xr - A")] 

where Xr is the actual value of the variable of the given system with the feedback control 
under consideration and A^ is the value of the coor dinate on the desired unstable 
limit cycle which is required to be stabilized. A" is obtained by allowing the system to 
evolve freely on the desired unstable orbit with the same parameter values but without 
imposing the control. 

We have implemented this idea for controlling unstable limit cycles of the Lorenz 
system. We choose parameter values u = 10, 6 = 8/3, r = 28. For these parameters one 
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of the unstable limit cycles has a point with coordinates (-12.786189, -19.364189, 24.00) 
[13]. The A of the system is found to be positive. Since the limit cycle is unstable, so 
starting with an initial state, even very close to the limit cycle, the trajectory will diverge 
away from the limit cycle. But after implementation of the control, multiplying parameter 
b by exp[e{Z — Z'^)], the Lyapunov exponents are found to be functions of e and for a 
certain range of values of e, the real part of all the Lyapunov exponents become negative 
thereby implying the stabilization of the limit cycle. Now, nearby trajectories are found to 
converge to the desired limit cycle as shown in fig. (3) also. The basin of attraction of the 
limit cycle is not infinite implying that there is a set of initial values starting from which 
the trajectories converge to the desired limit cycle. Otherwise they may escape to infinity. 
The transient r shows a similar trend as in the previous sections. But for a given accuracy 
CO the minimum value of r corresponding to the optimum e is much higher compared to the 
case of the unstable fixed point of the Lorenz system. Moreover, the variation of r with e 
for a given uj is not as rapid as in the case of the unstable fixed point. Consequently the 
T vs e plot is quite fiat and an optimum e cannot be obtained very accurately. 

We have also tested the control for another limit cycle of the Lorenz attractor. For the 
same parameter values, the coordinates of a point on a different limit cycle are (—13.917865, 
—21.919412, 24.00) [13]. The control was found to work effectively in this case with the 
same features as have been listed above. 

3.2 STABILIZATION OF HIGHER PERIOD ORBITS OF A 

DISCRETE DYNAMICAL SYSTEM 

We can extend the above control to convert unstable fixed points of higher 
period say {X*^, X*^, , X*^) to stable fixed points for given values of the control param- 
eters. What is required is a feedback that encodes as much information about the periodic 
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orbits as is necessary for its unique characterization. But there is a practical problem here. 
Since in this case the form of the control given by eq. (2) requires the convergence of the 
chosen variable to a set of k values of that variable, the controlling technique diminishes 
in utility with increase in period. For higher period orbits of a discrete dynamical system 
the more effective control is one which employs a logical OR structure in the feedback 
function. So in order to stabilize unstable fixed points of period k i.e. {X**} one of the 
parameters say fXj in (1) is multiplied by 

k 

exp[ell{{Xi)^-Xn] (10) 

i=l 

where {Xi)^ is the value of one of the chosen variables Xi at time = n of the modulated map 
after applying the control and (Xi*^ , Xi*^ , ....,Xi*'') is the set of k values of the variable 
Xi on the desired unstable period k orbit. 

We implemented this in the case of the logistic map for stabilizing period 2 orbit 
having X*^ and X*^ as the fixed points and found it to work effectively. 

4 STABILIZATION OF CHAOTIC ORBIT 

We have also tried to converge different trajectories to a particular desired 
chaotic trajectory using our control. Again taking the Lorenz system as an example we 
choose parameter values cr = 16, 6 = 4, r = 40. The real parts of the Lyapunov exponents 
are (1.37, 0.0, —22.37 ) [14] implying that this is a chaotic regime. We chose a chaotic 
trajectory starting with initial coordinates (10.0, 0.0, 30.0). We applied our control by 
multiplying b by exp[e{Z — Z'^)], where Z'^ is the Z coordinate of a point on the freely 
evolving chaotic trajectory. We found that the real parts of all the Lyapunov exponents 
of the modified Lorenz system become negative for a certain range of values of e, showing 
that the desired chaotic trajectory has become a stable trajectory and different closeby 
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trajectories starting from different initial states converge to the desired trajectory. For a 
given e the reciprocal of the slope of r vs ln(a;/a;o) is in good agreement with the minimum 
value of A. It has been verified numerically that the system settles down on to the desired 
trajectory in the presence of the control. The features of the control remain the same as 
in the previous section. We have also tested our control for other chaotic trajectories and 
found it to work effectively. 

5 STABILIZATION OF ARBITRARY FIXED POINTS FOR 
DISCRETE SYSTEM 

With exponential control it is possible to create new stable attractors which 
do not exist in the unmodulated system. This allows the modulated system to settle 
down to arbitrary fixed points which are not the fixed points (stable or unstable) of the 
unmodulated map. As an example we again take the logistic map given by eq. (4). Suppose 
the requirement is to stabilize the system to say a period 2 attractor (Xi, X2) which does 
not exist in the unmodulated map, implying that Xi and X2 are not the fixed points of 
F'^{li,X). The equations governing the dynamics in the presence of the control are 

= A^i^exp[e{Xr, - X,){Xn - X2)]X,(1 - X^) 

(11) 

Xn+2 = AlJ2exp[e{Xn+l - Xi){Xn+l - X2)]Xn+l{l - X^+l) 

where /ii and ii2 are given by 

//I =X2/(4Xi(l-Xi)) 

and 

fi2=Xi/ {4X2(1 -X2)) 

Any arbitrary combination of Xi and X2 cannot be stabilzed using eq. (11). Only those 
combinations of Xi and X2 can be stabilized in the logistic map using the exponential 
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control for which the A of the modulated map is negative. The coordinates of the points 
lying in the region enclosed by the curves and the coordinate axes in fig. (4) represent 
such combinations of Xi and X2 for which A is negative, and only such combinations can 
be stabilized using the map given by eq. (11). 

This procedure can easily be extended to stabilize the system to higher period orbits 
even when the points on the orbit do not correspond to either stable or unstable fixed 
points of the higher period orbits of the original system. 

6 CONCLUSION 

Exponential control is found to be effective both for maps as well as for flows. Since 
the knowledge of only a subset of variables on the desired unstable orbit is sufficient to 
settle the system on to that orbit, this makes the control more signiflcant particularly for 
systems with several degrees of freedom. In the presence of control, there exists a range 
of values of the stiffness constant e for which a repellor - an unstable fixed point, a limit 
cycle or a chaotic trajectory, of a system can be converted into an attractor. The A is 
found to be the function of e. The transient time r is a decreasing function of e for a given 
accuracy a; . However there exists an optimum stiffness control beyond which increasing e 
can increase r. This behaviour of r is found to hold for different types of orbits of different 
systems and is similar to the behaviour of A with e. The optimum e corresponds to the 
value of e where A is minimum. 

Exponential control is found to settle the Lorenz system on to its non-zero unstable 
fixed points faster than the control of Pecora and Carroll. In addition, it is also possible 
to stabilize arbitrary fixed points, chosen according to one's requirements, which do not 
correspond to either stable or unstable fixed points of the unmodulated map. 

We are investigating the effect of exponential control in the presence of noise. It 
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would be interesting to see whether such control is capable of controlling spatio-temporal 
chaos in coupled map lattice systems [15]. It would also be interesting to see whether 
exponential control can successfully recover the original signal when the original signal has 
been masked by the addition of noise [16]. 
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FIGURE CAPTIONS 



Fig.l The variation of A with e for the Lorenz system while stabihzing the unstable fixed 
point {■\/h{r — 1), ^/h{r^^V),r — 1) using the exponential control involving the parameter 
b and the variable Z, with cr = 10, 6 = 8/3 and r = 60. 

Fig.2 Plot of the transient time r vs ln(a;/a;o) for the Lorenz system with a = 10, b = 8/3 
and r = 60 while stabilizing the unstable fixed point {^/b{r — 1), \/b{r — 1), r — 1) using 
Pecora and Carroll's method with X = y^b{r — 1) drive (triangles), with Y = ^/b(r'^^l) 
drive (squares), and also with exponential control involving the parameter r and the vari- 
able X (circles), and r and Y (stars). 

Fig. 3 The dotted trajectory shows the convergence of a nearby trajectory to the unstable 
closed orbit, represented by the bold curve, projected on to the (X,Z) plane of the Lorenz 
system with cr = 10, 6 = 8/3, and r = 28 and a point with coordinates (—12.786189, 
-19.364189, 24.00). 

Fig. 4 The coordinates of the points lying in the region enclosed by the plotted curves and 
the coordinate axes represent the combinations of Xi and X2 which can be stabilized in 
the logistic map using exponential control for e = 0.1 and 1.0. 
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TABLE CAPTION 

The values of A, the largest real part of the Lyapunov exponents, of the Lorenz attractor 
under exponential control while stabilizing the fixed point {■\/h{r — 1), ■\/h{r — 1), r — 1) 
for cr = 10, 6 = 8/3 and r = 60 for different combinations of parameters and variables. 
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TABLE 1 





A for 


Control 
form 


e = -0.001 


e = -0.01 


e = -0.1 


r-X 


0.616 


-0.267 


-3.950 


r-Y 


0.394 


-2.381 


-7.180 


r-Z 


0.767 


1.117 


2.2.50 


b-X 


0.739 


0.927 


2.425 


b-Y 


0.707 


0.611 


-0.113 


b-Z 


0.786 


1.403 


7.884 


a-X 


0.717 


0.717 


0.717 


a-Y 


0.717 


0.717 


0.717 


a-Z 


0.717 


0.717 


0.717 
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